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METHOD FOR EXCLUSION OF ECTOPIC EVENTS FROM HEART RATE 

VARIABILITY METRICS 

Related Applications 

5 This application is a continuation-in-part of the following co-pending assigned 

patent applications: U.S. Patent Application Serial No. 09/704,844, filed November 2, 
2000; U.S. Patent Application Serial No. 09/802,316, filed March 8, 2001; U.S. Patent 
Application Serial No. 10/035,009, filed December 28, 2001; and U.S. Patent Application 
Serial No. 10/436,876, filed May 12, 2003, the disclosures of which are incorporated by 
10 reference in their entirety. 

Field of the Invention 

This invention pertains to cardiac rhythm management devices such as 
pacemakers and implantable monitoring devices. 

15 

Background 

Heart rate variability (HRV) refers to the changes in the length of time between 
consecutive heart beats during sinus rhythm and is primarily due to the interaction 
between the sympathetic and parasympathetic arms of the autonomic nervous system. 

20 Measurement and analysis of heart rate variability is thus a useful and non-invasive tool 
for assessing the status of the autonomic nervous system. 

The interval between successive heart beats, referred to herein as a BB interval, 
may be measured from a surface ECG or intracardiac electrogram as the time from the 
peak of one R wave representing ventricular depolarization to the peak of the next, 

25 referred to as an RR interval, or as the time from the peak of one P wave representing 
atrial depolarization to the peak of the next, referred to as an PP interval. Other measures 
of BB interval are possible and the present subject matter is not limited to the measures of 
BB interval based on morphological peaks described above. The variability of normal 
BB intervals (i.e., during sinus rhythm) can be determined and analyzed in several 

30 different ways in either the time domain or the frequency domain. Time domain 
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measurements involve the computation of statistics based upon the individual BB 
intervals making up a BB time series such as the standard deviation of the BB intervals in 
the series. Frequency domain analysis employs methods such as the Fast Fourier 
Transform (FFT) or autoregressive analysis to analyze the frequency spectrum of the 
5 variability in the BB intervals. Analysis of the frequency spectrum of heart rate 
variability has proven to be particularly valuable in assessing the relative activities of the 
sympathetic and parasympathetic nervous systems in a subject. Such assessment of the 
state of autonomic balance is a useful function for implantable cardiac rhythm 
management devices such as pacemakers and implantable cardioverter/defibrillators to 

10 perform as it can be used to modify the manner in which therapy is delivered by the 
device or to predict the occurrence of arrhythmias. 

Heart rate variability is a useful indicator of a subject's physiological state 
because it reflects the influence of internal feedback mechanisms within the body on the 
intrinsic heart rate produced at the sino-atrial node of the heart. Such feedback 

15 mechanisms are based on the interplay of the sympathetic and parasympathetic branches 
of the autonomic nervous system as well as various hormonal responses. Ectopic cardiac 
activity, such as premature atrial contractions (PAC's) and premature ventricular 
contractions (PVC's), are believed to be independent of the aforementioned internal 
feedback mechanisms. A heart rate variability metric calculated from a BB time series in 

20 which ectopic beats are present will be corrupted and may not accurately reflect the 
subject's true physiological state. For purposes of heart rate variability analysis, it is 
therefore desirable to exclude those BB intervals from a BB time series which are due to 
ectopic cardiac activity and are not the result of a normally conducted heart beat. It is 
with this problem that the present invention is primarily concerned. 

25 

Summary 

The present invention relates to an interval-based technique for removing ectopic 
or non-sinus events from a BB interval signal which may be implemented in an 
appropriately programmed implantable medical device. The technique involves filtering 
30 the BB intervals to exclude those BB intervals which are greater or less than a function of 
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preceding BB intervals by a specified threshold value, the excluded BB intervals being 
deemed to be due to ectopic cardiac activity. Such a function of preceding BB intervals 
may be a statistic (e.g., a median or an average) computed from a plurality of preceding 
BB intervals. A heart rate variability metric may then be calculated from the filtered BB 
5 intervals. In one embodiment, the filtering is performed by comparing the present BB 
interval to a statistic computed from the plurality of previous BB intervals stored in a 
first-in-first-out buffer, where the buffer contains a maximum number N of preceding BB 
intervals. The buffer is then updated by removing the oldest interval and storing the 
present BB interval therein if the present BB interval was not excluded as ectopic. 

10 

Brief Description of the Drawings 

Fig. 1 is a block diagram of an exemplary cardiac rhythm management device for 
practicing the present invention. 

Fig. 2 illustrates an exemplary algorithm for filtering BB intervals to remove 
15 ectopic beats. 

Fig. 3 is a graph showing the results of interval-based ectopic filtering of a 
representative RR interval series. 

Detailed Description 

20 As noted above, heart rate variability refers to the variability of the time intervals 

between successive heart beats during a sinus rhythm. Spectral analysis of heart rate 
variability involves decomposing a signal representing successive beat-to-beat intervals 
into separate components representing the amplitude of the signal at different oscillation 
frequencies. It has been found that the amount of signal power in a low frequency (LF) 

25 band ranging from 0.04 to 0.15 Hz is influenced by the levels of activity of both the 
sympathetic and parasympathetic nervous systems, while the amount of signal power in a 
high frequency band (HF) ranging from 0.15 to 0.40 Hz is primarily a function of 
parasympathetic activity. The ratio of the signal powers, designated as the LF/HF ratio, is 
thus a good indicator of the state of autonomic balance, with a high LF/HF ratio 

30 indicating increased sympathetic activity. If an implantable medical device is 
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programmed to monitor the LF/HF ratio, the device can log a clinically significant event 
when the ratio exceeds a specified threshold value, as well as possibly automatically 
altering its mode of operation (e.g., initiating different therapies or performing more 
computationally intensive data analysis for arrhythmia prediction). 
5 A series of BB interval values can be regarded as a discrete signal indexed by 

heartbeat such that each value of the signal represents a BB interval for a particular 
heartbeat. In order to properly analyze the frequency content of heart rate variability, 
however, the BB time series can be re-sampled at a specified sampling frequency so as to 
equalize the time intervals between interval values and thus convert the time series into a 

10 discrete time signal, where the sampling frequency is selected to meet the Nyquist 
criterion with respect to the frequencies of interest. Spectral analysis of such a BB 
interval signal can then be performed directly in the frequency domain using discrete 
Fourier transform or autoregression techniques in order to compute the LF/HF ratio. A 
time-domain technique for determining the high and low frequency components of the 

15 signal could also be used in which the BB interval signal is input to low band and high 
band digital filters and signals proportional to the power of the BB interval signal in each 
of the low frequency and high frequency bands are derived so that the LF/HF ratio may be 
computed. 

Both frequency domain and time domain analysis performed as described above 
20 are computationally intensive, however, and require the storage of large amounts of BB 
interval data. Such methods may therefore not be practical in a typical implantable 
medical device which is a small battery-powered device with limited processing power. 
Statistical techniques for heart rate variability analysis that do not involve such processing 
overhead are described in above-cited U.S. Patent Application Ser. No. 10/436,876, 
25 where it should be appreciated that the term RR interval in that application can be taken 
to mean any measure of the interval between heart beats, referred to herein as a BB 
interval. Such techniques are used to generate surrogate parameters which represent 
different frequency components of a BB time series and from which the LF/HF ratio or 
other quantity may be calculated. The rMSSD statistic, for example, is defined as the 
30 square root of the mean of the squared successive differences of a BB time series: 
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rMSSD = E{(BBi-BBj_i) 2 } 05 

where E is the expectation or mean value operator, and BBj refers to the i th BB interval in 
5 the series. The square root step in the calculation can be omitted to give the [rMSSD] 2 
parameter. By averaging the square of the successive interval-to-interval difference 
values in the BB time series, the rMSSD or [rMSSD] 2 statistic maximally reflects 
variations in the BB intervals that occur with each successive interval and progressively 
attenuates variations in the BB intervals that occur at lower frequencies. The frequency 

10 response represented by the rMSSD statistic is therefore greatest at the maximum 
frequency that can be represented in the time series (i.e., approximately one-half of the 
average heart rate, similar to the Nyquist frequency in a time series with regular intervals) 
and then decreases linearly with decreasing frequency so that lower frequency variability 
in the BB time series is not represented. Computation of the rMSSD or [rMSSD] 2 

15 statistic thus captures a frequency range of heart rate variability which is similar to the HF 
band. The SD 7 statistic is defined as the standard deviation of the mean values of all 
successive 7-second segments in the BB time series, and [SD7] 2 is the square of that 
standard deviation or variance. By computing mean values of the BB intervals over 7- 
second segments, the [SD 7 ] 2 statistic averages out the variations in the BB intervals that 

20 occur over time intervals equal to or smaller than 7 seconds (i.e., variations at frequencies 
higher than approximately 0.15 Hz), with the variations in the BB intervals occurring 
over longer intervals then being reflected by computation of the variance of those mean 
values. The [SD 7 ] 2 statistic thus captures those frequency components of the BB interval 
signal from some frequency above DC (since a variance calculation eliminates the DC 

25 component of a signal) to approximately 1/7 or 0.15 Hz (i.e., a frequency with a period 
equal to the length of the 7-second segment). The SD25 statistic is similarly defined as 
the standard deviation of the mean values of all successive 25-second segments in the BB 
time series, where [SD25] 2 is the square of that standard deviation or variance. The 
[SD25] 2 statistic thus captures those frequency components of the BB interval signal from 

30 some frequency above DC to approximately 1/25 or 0.04 Hz (i.e., a frequency with a 
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period equal to the length of the 25-second segment). If the [SD25] 2 statistic is subtracted 
from the [SD 7 ] 2 statistic, the resulting parameter captures a frequency range of heart rate 
variability which is similar to the LF band. 

The rMSSD, [SD25] 2 , and [SD 7 ] 2 statistics may thus serve as surrogates for the 
5 frequency components of a BB time series. An estimate of the LF/HF ratio may then be 
computed as: 

Estimated LF/HF = K {[SD 7 ] 2 - [SD 25 ] 2 }/ [rMSSD] 2 

10 where K is a constant. A linear regression analysis may be performed in which the 
estimated values are correlated with the actual spectrum of a BB time series to derive the 
value of K. Alternatively, the estimated LF/HF ratio may be compared with appropriately 
scaled threshold values in order to assess the autonomic balance of a subject, which 
eliminates the need for K in the calculation. 

15 Heart rate variability metrics thus include parameters computed by time-domain 

filtering of a BB interval signal, parameters computed by frequency-domain analysis of an 
BB interval signal, and statistical parameters which serve as surrogates for different 
frequency components of the BB interval signal. Regardless of the particular heart rate 
variability metric, however, it is desirable to exclude from the BB interval signal those 

20 BB intervals which are due to ectopic cardiac activity such as PVC's. The present 
invention relates to a BB interval filtering technique for removing non-sinus rhythm 
events from a BB interval series used for heart rate variability analysis. A description of 
the technique is set forth below following a description of an exemplary hardware 
platform. 

25 

1 . Exemplary implantable device description 

The present invention may be incorporated into any cardiac device with the 
capability of sensing cardiac electrical activity, including devices for monitoring only and 
those for delivering therapy in the form of electrical stimulation to the heart. For 
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illustrative purposes, however, the invention will be described with reference to a dual- 
chamber pacemaker (i.e., one that senses and paces both the atria and ventricles). 

Cardiac rhythm management devices are implantable devices that provide 
electrical stimulation to selected chambers of the heart in order to treat disorders of 
5 cardiac rhythm. A pacemaker, for example, is a cardiac rhythm management device that 
paces the heart with timed pacing pulses. The most common condition for which 
pacemakers are used is in the treatment of bradycardia, where the ventricular rate is too 
slow. Cardiac rhythm management devices may also treat tachyarrhythmias, where the 
heart rate is too fast, by anti-tachycardia pacing and/or delivery of defibrillation shocks. 

10 Such devices are usually implanted subcutaneously on the patient's chest and connected 
to electrodes by leads threaded through the vessels of the upper venous system into the 
heart. An electrode can be incorporated into a sensing channel that generates an 
electrogram signal representing cardiac electrical activity at the electrode site and/or 
incorporated into a pacing or shocking channel for delivering pacing or shock pulses to 

15 the site. 

A block diagram of an implantable cardiac rhythm management device is shown 
in Fig. 1. The controller of the device is made up of a microprocessor 10 communicating 
with a memory 12 via a bidirectional data bus, where the memory 12 typically comprises 
a ROM (read-only memory) and a RAM (random-access memory). The controller could 

20 be implemented by other types of logic circuitry (e.g., discrete components or 
programmable logic arrays) using a state machine type of design, but a microprocessor- 
based system is preferable. As used herein, the programming of a controller should be 
taken to refer to either discrete logic circuitry configured to perform particular functions 
or to executable code stored in memory or other storage medium. The controller is 

25 capable of operating the device so as to deliver a number of different therapies in 
response to detected cardiac activity. A telemetry interface 80 is also provided for 
enabling the controller to communicate with an external programmer. 

The embodiment shown in Fig. 1 has two sensing/pacing channels, where a 
pacing channel is made up of a pulse generator connected to an electrode while a sensing 

30 channel is made up of the sense amplifier connected to an electrode. A MOS switch 
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matrix 70 controlled by the microprocessor is used to switch the electrodes from the input 
of a sense amplifier to the output of a pulse generator. The switch matrix 70 also allows 
the sensing and pacing channels to be configured by the controller with different 
combinations of the available electrodes. The channels may be configured as either atrial 
5 or ventricular channels. In an example configuration, an atrial sensing/pacing channel 
includes ring electrode 43a and tip electrode 43b of bipolar lead 43c, sense amplifier 41, 
pulse generator 42, and a channel interface 40. A ventricular sensing/pacing channel 
includes ring electrode 33a and tip electrode 33b of bipolar lead 33c, sense amplifier 31, 
pulse generator 32, and a channel interface 30. The channel interfaces communicate bi- 

10 directionally with a port of microprocessor 10 and may include analog-to-digital 
converters for digitizing sensing signal inputs from the sensing amplifiers, registers that 
can be written to for adjusting the gain and threshold values of the sensing amplifiers, and 
registers for controlling the output of pacing pulses and/or changing the pacing pulse 
amplitude. A shock pulse generator 20 is also interfaced to the controller for delivering 

15 defibrillation shocks through electrodes selected by the switch matrix. In the illustrated 
embodiment, the device is equipped with bipolar leads that include two electrodes which 
are used for outputting a pacing pulse and/or sensing intrinsic activity. Other 
embodiments may employ unipolar leads with single electrodes for sensing and pacing. 
The switch matrix 70 may configure a channel for unipolar sensing or pacing by 

20 referencing an electrode of a unipolar or bipolar lead with the device housing or can 60. 

The controller 10 controls the overall operation of the device in accordance with 
programmed instructions stored in memory. The controller 10 interprets electrogram 
signals from the sensing channels in order to control the delivery of paces in accordance 
with a pacing mode and/or deliver shock therapy in response to detection of a 

25 tachyarrhythmia such as ventricular fibrillation. The sensing circuitry of the device 
generates atrial and ventricular electrogram signals from the voltages sensed by the 
electrodes of a particular channel. An electrogram is analogous to a surface ECG and 
indicates the time course and amplitude of cardiac depolarization that occurs during either 
an intrinsic or paced beat. When an electrogram signal in an atrial or sensing channel 

30 exceeds a specified threshold, the controller detects a chamber sense (i.e., an atrial or 
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ventricular sense, respectively), which pacing algorithms may employ to trigger or inhibit 
pacing and from which heart rates may be derived by measuring the intervals between 
senses. The measured intervals between successive ventricular senses or between 
successive atrial senses make up a BB time series. As described below, the controller 
5 may also be programmed to eliminate ectopic beats from the BB time series used to 
calculate heart rate variability metrics. 

2. Interval-based filtering of ectopic events 

A cardiac rhythm management device such as illustrated in Fig. 1 can be 

10 programmed to determine heart rate variability by analyzing data received from its 
sensing channels. The intervals between successive cardiac chamber senses, referred to 
as BB intervals, can be measured for a specified period of time or a specified number of 
beats and their variability analyzed. A typical BB time series, for example, would be 
made up of BB intervals over 24 hours or other period. As described above, a BB 

15 interval signal collected in this manner can be spectrally analyzed by the device in order 
to determine the frequency content in the LF and HF bands by either transforming the 
signal into the frequency domain, decomposing the signal with bandpass filters, or by 
calculating statistical parameters which serve as surrogates for the actual specific 
frequency components. 

20 One way to derive a BB interval signal representing heart rate variability during a 

sinus rhythm is to remove ectopic beats with an event-based technique. For example, 
ectopic ventricular beats (i.e., premature ventricular contractions or PVCs) or ectopic 
atrial beats (i.e., premature atrial contractions or PACs) can be detected by monitoring 
whether an atrial sense (or P wave) precedes each ventricular sense (R wave). If two 

25 ventricular sense events occur with no atrial sense event between, the second ventricular 
sense is designated as a PVC. BB intervals adjacent to the PVC can then be removed 
from the BB interval series used for heart rate variability calculations (or changed to 
interpolated values). A similar approach can be applied when consecutive atrial senses 
are detected with no ventricular sense between. There are limitations with event-based 

30 ectopic beat removal, however. For example, ectopic events may be missed when a PVC 
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following an atrial sense is deemed a sinus rhythm beat or when a PAC that is conducted 
through to the ventricle is classified as a sinus rhythm beat. Also, it is difficult to 
evaluate the performance of the implantable device based upon surface ECG Holter data 
because of the limited visibility of atrial events. In accordance with the present invention, 

5 a BB interval-based ectopic removal algorithm is used either instead of or in addition to 
the event-based technique described above. With this technique, BB intervals are 
classified as ectopic or not based on a history of previous BB intervals. An implantable 
medical device may thus be programmed to perform heart rate variability analysis by: 1) 
sensing electrical activity in a cardiac chamber and generating a chamber sense signal 

10 when the sensed electrical activity exceeds a predetermined threshold; 2) measuring time 
intervals between each pair of successive chamber senses, referred to as BB intervals, the 
measured BB intervals constituting an BB time series; 3) filtering the BB intervals to 
exclude those BB intervals which are greater or less than a statistic computed from a 
plurality of preceding BB intervals by a specified threshold value, the excluded BB 

15 intervals being deemed to be ectopic BB intervals which are due to ectopic cardiac 
activity; and, 4) computing a heart rate variability metric from the filtered BB intervals. 
The statistic computed from the plurality of BB intervals may be, for example, a median 
or a weighted average, and the threshold value may be a specified percentage of the 
calculated statistic. The device may also be programmed to maintain a count of ectopic 

20 events detected by the interval-based filtering or in conjunction with an event-based 
technique for diagnostic purposes. 

Fig. 2 illustrates an exemplary algorithm for removing ectopic events in 
accordance with the invention. Successive non-ectopic BB intervals are stored in a first- 
in-first-out (FIFO) buffer, and a median value of the stored intervals is computed for 

25 evaluating present BB intervals as they are collected. The buffer has a maximum capacity 
of N intervals where N is a specified integer. In a presently preferred embodiment, for 
reasons of both performance and computational efficiency, the number N is specified as 
three. Starting at step 201, the device acquires the next present BB interval. At step 202, 
the present BB interval is compared with a limit value HighLim and a low limit value 

30 LowLim. If the BB interval is above the high limit value or below the low limit value, 
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the interval is classified as ectopic. Exemplary values for the high and low limit values 
are 2000 ms and 300 ms, respectively. Thus, if the present BB interval is not within the 
allowed range, step 206a classifies present BB interval as ECTOPIC_OUT_OF_RANGE, 
and the interval is not used for HRV metric computation. The oldest BB interval in the 
5 buffer is also removed, and the device returns to step 201 to collect the next BB interval. 
If the present BB interval is within the specified higher and lower limit values, the buffer 
is checked at step 203 to see if it contains N intervals. If so, the median MD of the N 
intervals in the buffer is computed at step 204b, and the present BB interval is compared 
with the computed median at step 205. If the present BB interval differs from the median 

10 MD by more than a specified threshold value Th (which may be a percentage of the 
computed median), the algorithm proceeds to step 206a so that the interval is classified as 
ECTOPIC J3UTJ3F_RANGE. If not, the interval is classified as OK at step 206b so 
that it can be used in the calculation of a heart rate variability metric. Next, the buffer is 
updated at step 207 by: 1) if the buffer contains N intervals, removing the oldest BB 

15 interval from buffer, and 2) storing the present BB interval in the buffer. The algorithm 
then returns to step 201. Referring back to step 203, if the buffer does not contain N 
intervals, the present BB interval is not compared with a median computed from the 
intervals contained in the buffer. Instead, at step 204a, the present BB interval is 
classified ECTOPIC JNSUFFICffiNTJNTERVALS (i.e., indeterminate) and not used 

20 for HRV metric computation. The algorithm then updates the buffer at step 207. 

Fig. 3 is a graph showing the results of performing interval-based ectopic filtering 
of ectopic events as described above on a representative RR interval series. In this 
example, the statistic computed from preceding RR intervals is the median of the 
previous three non-ectopic RR intervals. An RR interval is classified as ectopic if it is 

25 greater or less than the median by a specified percentage of the computed median, shown 
graphically as those RR intervals which are above a maximum value Max or below a 
minimum value Min. Two such RR intervals are shown in the figure as being classified 
as ectopic. 
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Although the invention has been described in conjunction with the foregoing 
specific embodiments, many alternatives, variations, and modifications will be apparent 
to those of ordinary skill in the art. Other such alternatives, variations, and modifications 
are intended to fall within the scope of the following appended claims. 
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